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BLOCK-DIAGONALIZATION OF ODES IN THE SEMICLASSICAL LIMIT 
AND C u VS. C°° STATIONARY PHASE 

OLIVIER LAFITTE, MARK WILLIAMS, AND KEVIN ZUMBRUN 


Abstract. Motivated by issues in detonation stability, we study existence of block-diagonalizing 
transformations for ordinary differential semiclassical limit problems arising in the study of high- 
frequency eigenvalue problems. Our main results are to (i) establish existence of block-diagonalizing 
transformations in a neighborhood of infinity for analytic-coefficient ODE, and (ii) establish by a 
series of counterexample sharpness of hypotheses and conclusions on existence of block-diagonalizing 
transformations near a Enite point. In particular, we show that, in general, bounded transformations 
exist only locally, answering a question posed by Wasow in the 1980’s, and, under the minimal 
condition of spectral separation, for ODE with analytic rather than C°° coefficients. The latter 
issue is connected with quantitative comparisons of C" vs. C°° stationary phase estimates. 


1. Introduction 

Motivated by problems in detonation and related hydrodynamical and continuum-mechanical 
stability, we consider the general semi-classical limit problem 

(1.1) h(d/dx)Z = (A(x,h;q) + h,B(x,h;q))Z, Z € C N , h . —>• 0 + , 

on a possibly unbounded domain x € [a, b] C M, representing a generalized spectral problem with 
wavelength h € R + and frequency k = 1/h. Here, q € R s , bounded, records any additional 
parameters associated with the problem: typically, spectral angle and or bifurcation parameters. 
Such systems arise for example as generalized eigenvalue problems for the linearized equations 

(1.2) m t = Lv:=Y^ a^x^v, d “ := d" * 1 • • • d^, 

a 

about a steady planar solution 

(1.3) w(xi,...,x d ,t)=w(xi) 

of a PDE 

(1.4) w t = T(w) := X f a (w)d£w, 

a 

in the high-frequency limit. Specifically, taking the Fourier transform in (x 2 , ■ ■ ■ ,x d ), reduces the 
eigenvalue problem Xv = Lv to an ODE 

(1.5) Xv = L^v := 53o Q (xi)^i 1 (i6)“ 2 ••• (i£d) ad - 

a 

Writing (11.51) as a first-order system in a suitable phase variable Z including v and appropriate 
additional xi-derivatives, and defining h~ l := |£,A| as spectral frequency and 

? = (a ,i) ■■= (e/|A,ei,e/|A,ci) 
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as spectral angle, we arrive in the high-frequency limit |£, A| —>• oo at a problem of form (11.11) 
in the variable x = x±. In this context, values (q,h) for which there exist solutions Z of (11.ID 
satisfying appropriate boundary conditions at endpoints x = a, b correspond to spectra A = X/h 
of the Fourier transform Lg of the linearized operator about the wave, hence an understanding of 
small- h behavior of ED corresponds to an understanding of high-frequency spectral stability. See 
[Erl IZ11ILWZ11ILWZ21 for specific examples pertaining to stability of detonation waves. 

Our goal in this paper is a systematic treatment of local block reduction of ED, or decomposition 
of the equations into spectrally separated blocks possessing nontrivial turning points, in particular 
in the important case, not previously treated to our knowledge, of a neighborhood of plus or minus 
infinity. At finite points, for which existence of locally diagonalizing transformations has been 
exhaustively studied in EOT, our goal is to determine sharpness of hypotheses and conclusions, 
and in particular compare results obtainable by complex-analytic methods to those obtained by C r 
methods in, e.g., [ ZlilLWZlj . The treatment of the resulting smaller blocks after this decomposition, 
and the global implications for stability are studied, for example, in [W1 fOl LWZ2 , 1LWZ3II . 

For simplicity of exposition, we will suppress in what follows dependence on the parameter q 
(corresponding in example (11.51) to restriction to the ID case d = 1), leaving only dependence of 
coefficients on the parameter h. However, it is an important point that all estimates of the paper 
carry over to the general case, uniformly in the parameters q and h, the treatment of g-dependence 
being no different than the treatment of dependence on h. As noted in [LWZlj . such uniform 
estimates are important in the verification of high-frequency stability, or nonexistence of unstable 
spectra for h sufficiently small, a property involving all parameter values, as opposed to instability, 
a property that need be checked only at isolated strategically chosen parameter values. 

This, and the treatment for unbounded domains of block-diagonalization at infinity, are two of the 
main goals of the present analysis. In the companion paper |LWZ2] , we have already made good use 
of these ideas, applying and further extending then£| to obtain the result of high-frequency stability 
of detonation waves in certain media, making rigorous the important observations of Erpenbeck [Er] 
made by a combination of formal and rigorous analysis in the 1960’s, but up to now not rigorously 
verified. At the same time, given the sometimes bewildering array of different techniques that 
have been developed for the analysis of this problem, including asympotic ODE and microlocal 
analysis/WKB expansion, and C r vs. analytic stationary phase, we seek to make clear what can 
and cannot be accomplished under various assumptions on (II.ID : that is, to remove the uncertainty 
whether a stronger result could perhaps be obtained by a different technique. 

Here, our main result is to make an explicit connection between existence of block-diagonalizing 
transformations and decay rates for certain oscillatory integrals 

(1.6) lim / e T > a(y,h)dy, 

h-+ o+ J a 

whereby we are able to resolve a number of such questions by stationary phase computations under 
appropriate conditions on the symbol a. In particular, we show that: (i) block-diagonalization can 
in general be done only locally, answering a longstanding question posed by Wasow in his 1985 text 
[W] , and (ii) in general requires analyticity and not just C r or C°° of the coefficients of ED- 

The former is discussed in the Remark, p. 89 of jWj . comparing analogs of our Theorem ll.3l to an 
analog of our Theorem ll.il “Theorems 6.1-1 and 6.1-2 are strictly local, although the decomposition 
described in Theorem 12.3-1 of the Appendix is globally valid in all of D. A global uncoupling of 
the given differential equation by one and the same transformation with an asymptotic series in 
powers of epsilon, valid in large regions would be a boon to the theory. On the other hand it is quite 
possible that such a theorem does not exist, and then one would like to see counterexamples. ” 


1 Among other things, analyzing nontrivial turning points, finite and infinite, the case to which we here reduce. 
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The latter appears to be linked to interesting recently observed phenomena in spectral theory 
[HS] (almost-sure diffusion of spectra under random C°° perturbation of an analytic-coefficient 
operator) and propagation of singularities [Leb] (diffraction by C°° vs. analytic boundary in M 3 ). 
It is obtained via sharp stationary phase estimates (11.171) (11.181) for Gevrey class symbols a, inter¬ 
polating between the algebraic van der Korput bounds for C r symbols and the exponential bounds 
for analytic a; itself of independent interest, this estimate too so far as we know is new. 

Notation. Symbols ~, <, > indicate equality/inequality up to a constant factor bounded uni¬ 
formly with respect to parameters. a(M) indicates spectrum of a matrix or linear operator M. 

1.1. Background/previous results. Before stating our main results, we set the stage with a 
brief further discussion of some background and motivation for the analysis. 


1.1.1. Block diagonalization and WKB expansion. The classical WKB approach to approximating 
solutions of a system 

(1.7) hZ' = A(x)Z + hB(x)Z 


falling under the general form m is to seek a basis of approximate solutions of form 
.8) Zj(x) = e^M+ELo h*ki(x)ph,m( x ^ ph,m {x) = p. o(x) + hPj ^ x ) + . . . + h m P,,m(x), 

(hd x — <&(x, h))Zj = 0(h m \Zj\)- equivalently, a parametrix = p h e ( hH + K )(y )IS , p = diag {hj}, 
K = diag {Yli=o h- 1 where aj = d x hj and Pj$ are eigenvalues and eigenvectors of A(x), and 

ph,m ._ _ (Pn’ m ) denotes the matrix with columns p^ m . So long as the eigenvalues aj 

remain distinct, so that P h , ( P h ) -1 may be taken uniformly bounded, one can convert the formal 
h m modeling error to a rigorous convergence bound by a Lyapunov-Perron type integral equation 
mimicking the usual construction of invariant manifolds, in which jth parts of the propagator 
are integrated against modeling error along “progressive contours” for which 9f l(hj — h m )(z) is 
nonincreasing for all m jwi mi- These may be real contours if the eigenvalues of A maintain a 
neutral pairwise spectral gap, but in general are complex, requiring A. B analytic. 

Difficulties occur at nontrivial turning points, where eigenvalues of A collide, and, for unbounded 
domains, at oo, where the usual prescription of progressive contours breaks down. Here, we will seek 
not to carry out a complete expansion as in (11.81) . but only a block-diagonalization corresponding 
to invariant subspaces with distinct spectra of A: essentially a vector version of (11.81) . This is of 
course also a preliminary step to full conjugation, decoupling the problem into scalar modes and 
irreducible m x m blocks, m > 1, containing nontrivial turning points, to be analyzed by more 
special techniques as in, e.g., E2HDJ- Importantly, this includes block-diagonalization at infinity. 


1.1.2. C r vs C u diagonalization, and diagonalization on unbounded domains. A new aspect of the 
1- and multi-D analyses of detonation stability in [Zll ILWZlj was to carry out rigorous WKB-type 
expansion for (ElD on unbounded domains. This involved (i) effectively resumming the usual series 
expansion to obtain an exact solution at infinity/integrability of modeling error in x, and (ii) the 
development in |Zl| of a new “variable coefficient gap lemma” for C r -coefficient ODE, generalizing 
to the semiclassical setting the standard gap lemma of |GZj . by which one may obtain solutions 
with desired behavior at infinity so long as the associated eigenvalues of A(x) are (a) semisimple 
and (b) satisfy a neutral numerical range condition roughly corresponding to a neutral spectral 
gap (separation of real parts) from other eigenvalues, with A and B converging to their limits at 
infinity at L 1 rate. The description is valid globally on the (possibly unbounded) interval where 
(a)-(b) hold. However, both conditions (a)-(b) fail in general for the class of problems arising in 
high-frequency stability of detonation waves, as studied, e.g., in (Erl ILWZ2] . 

This raises the questions (partially alluded to above) whether: (1) the conditions (a)-(b) are 
indeed necessary for the C' r -coefficient problem, or whether there could be provided by different 
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techniques a more general C r -coefficient theorem requiring only separation and not spectral gap 
of the eigenvalues of A, the “natural” condition needed for formal WKB expansion, and (2) the 
classical local diagonalization results of Wasow [W] (for the analytic-coefficient problem) are sharp, 
or whether they could be extended to a global result valid on the whole interval on which the 
eigenvalues of A(x) remain separated. It is these two questions, and the physical issues originating 
in detonation that prompt them, that are the primary practical motivations for our analysis. 

1.2. Main results. 

1.2.1. The profile problem, assumptions, and approximate block-diagonalization. Returning to the 
PDE problem fjl.3D - fjl.4jl . consider the commonly-occurring case of a front- or pulse-type solution 
lim^-i-oo w(z) = w±. Writing the standing-wave ODE as a first-order system 

(1.9) Z' = F(Z), 

in a phase variable Z consisting of w and appropriate derivatives, we find, so long as w± are 
nondegenerate hyperbolic equilibria, i.e., the Jacobians dF(w±) possess no center subspace, that 
the profile w consists of the projection onto the ui-coordinate of a profile Z, lim^ioo Z(z) = Z±, 
of (11.91) . which in turn corresponds generically to a transversal intersection of the stable manifold 
at Z + and the unstable manifold at Z_ of m- Assuming that the coefficients of (11.41) are C r , 
we obtain from this construction a profile w that is C r+l in x, and (by standard stable/unstable 
manifold theorems) converges exponentially in r derivatives to its limits w± as x —>• ±oo. Thus, 
for smooth coefficients f a € C r , C r smoothness and exponential convergence at ±oo are natural 
conditions to impose on A and B in ED- 

For analytic coefficients f a € C w , we obtain by the same construction a solution vj that is 
analytic for all 2 € R; however, already at the level of profiles, the situation is slightly more subtle 
at z —>• ± 00 . Namely, as we show in Theorem I A. II of Appendix lA.21 the stable manifold construction 
in the analytic coefficient case yields the much stronger result of existence/analyticity in a wedge 

■Hz > 0, |9z| < 

with exponential decay |tc(z)| < C(fj)e~^ z \, u,fj > 0, and similarly for the unstable manifold at 
z —>• — 00 . Analyticity of A and B on a strip around the real axis and wedges around plus and minus 
infinity, with exponential convergence as Hz —>• ± 00 , are thus natural assumptions for ED in the 
analytic coefficient case. This observation so far as we know is new; moreover, the strengthened 
assumptions at ±00 turn out to be essential for our treatment of exact block-diagonalization. 

Our first main result (following) is that, under the above natural assumptions, and assuming 
spectral separation of eigenvalues of the limiting coefficient matrix, there exist global approximately 
block-diagonalizing transformations to all orders, preserving the original assumptions. 

Theorem 1.1 (Global approximate diagonalization). Consider a general semiclassical limit prob¬ 
lem hW 1 = A(x,h)W, x € [a, b] C R, A e C r (x,h), such that the eigenvalues of Af, 0) may be 
divided into two groups separated uniformly in x, h for all x € R. Then, for each 1 < k < r, 
there exists a uniformly invertible change of coordinates W = T k (x,h)Wk, T € C r , yielding the 
approximately diagonalized system 

(1-10) 4) 1U+ !,*(“ ;i)l w, 

A k ,9j € C r ~ k , with F s+l = T s + 0(h s+1 ). If a = —00 or b = +00 and A is exponentially 
converging as x —>• ±00 in up to r derivatives, uniformly in x, h, then A k is exponentially converging 
and 6 k exponentially decaying in r — k derivatives as x —>• ± 00 , uniformly in h. If, moreover, A is 
analytic on a strip around the real axis and wedges around plus and minus infinity, with exponential 
convergence as (Rir — > ± 00 , then T k , A k , 9 k have these properties as well. 
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Proof. This is an immediate consequence of Proposition 12.11 and Remark 12.31 below. 


□ 


Remark 1.2. For cr(dF(z±)) real, the proof of Theorem IA.1I yields analyticity of profiles on half¬ 
planes Rz > M, Rz < — M, M » 1. This was shown in the scalar case in |LWZ21 Proposition 
4.1] using a direct, Implicit Function Theorem argument. The example Z € C, F(Z) = Z 2 — Z, 
Z(x ) = — tanh(x) shows that this result is sharp, as tanh(ir) = tan(r) has poles at r = 7r/2 + 2nj, 
j € Z. Theorem I A. II extends this to the system case, allowing more general applications; see for 
example Remark 2.1 |LWZ2l p. 9] on treatment of detonations with multi-component reactions. 

1.2.2. Local exact block-diagonalization. Complementing the global approximate diagonalization 
result of Proposition UTTl we have the following local, analytic exact diagonalization result, esentially 
a finite h-regularity, finite-accuracy version of Wasow’s theorem |Wj in the h-analytic case. 

Theorem 1.3 (Exact local diagonalization at a finite point). Given an ODE 

(1.11) h\V' = 2 22 ) w + h p QW, p> 1, 

An, © uniformly analytic in x in a complex neighborhood of x = x* and continuous in h in a postive 
real neighborhood of h = 0, with no eigenvalues of Ajj in common, there exists a coordinate change 

W = T(x. h)Z. T — W ’“ 12 ) , such that hZ' = (' 411 + f 0u ^ Z, with a „ 

/3j uniformly analytic in x and continuous in h in neighborhoods of h = 0, x = x *, 

Proof. See Section l2~2l □ 

An important extension of Theorem 11.31 for the applications we have in mind is the following 
result giving existence of block-diagonalizing conjugators near infinity for a linear system (11.111) 
with Aj, 0 are analytic in x on the wedge Wm,/3 ■ > 0, |9x| < /3Rx, with 

(1.12) \Aj(x,h) -Aj(+oo,h)\ < Ce~^ x , |0(x, h)\ < C on Wp, 
for some limiting value Aj(+oo) and constant C, uniformy in h « 1. 

Theorem 1.4 (Exact block diagonalization at infinity). For Ajj, © as in (11.121) . if limits Ajj(+oo) 
have no eigenvalues in common, then there exists a uniformly bounded analytic conjugator T(-) 
on possibly smaller wedge Wm',P', for ff sufficiently small and M' > 0 sufficiently large, with 
| T — Id| = 0(h p ). Moreover, if 0 is uniformly exponentially convergent as Re{x) —>• oo, then T is 
as well, and if 0 is uniformly exponentially decaying as 9ft(x) —>• oo, then (T — Id) is also. 

Proof. See Section [2731 □ 

So far as we know, both statement and proof of Theorem 11.41 are new; indeed, as mentioned 
earlier, the novel assumption of analyticity on a wedge appears to be essential. A related problem 
is existence of block-diagonalizing conjugators for a singular ODE system 

(1.13) zh^- = ( A ^ .° \ (z, h)W + h p O(z, h)W 

dz \ U A 22 J 

in the vicinity of the singular point z = 0. Such problems arise, for example, in the treatment 
of “hybrid resonance” or “X-mode” heating of fusion plasma |DIWl iDILj . where coincidence of 
regular-singular and turning points lead to interesting physical phenomena. A first step in their 
rigorous analysis is block reduction to a 2 x 2 system equivalent to a modified Bessel equation 0. 

Corollary 1.5 (Exact block diagonalization at a singular point). Let Ajj, 0 in ((1.13D be analytic 
in z and continuous in h on B(0,r) x (0, ho) C C x R + and 4n(0), ^ 22 ( 0 ) have no eigenvalues 
in common. Then, there exists a uniformly bounded block-diagonalizing conjugator with uniformly 
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bounded inverse T = Id + 0{hP) of (11.131) . analytic in z on a slit ball around z = 0 with branch cut 
along the negative real axis, and continuous at z = 0. 

Proof. This follows from Theorem 11.41 via the transformation z —>• x = — In z\ see Section 12.51 □ 

Remark 1.6. Example 14.31 below shows that the conclusions of Corollary 11.51 are sharp even for Ajj, 
0 independent of h and analytic at z = 0; specifically, there need not then exist a diagonalizer that 
is analytic on a neighborhood of z = 0, despite a formal power series solution to all orders. At 
the same time, the hypotheses of analyticity in z on £>(0, r) may be weakened to analyticity in a 
neighborhood of the origin on the Riemann surface with branch cut along the negative real axis. 


1.2.3. Oscillatory integrals and counterexamples. Consider the 2x2 triangular system 

(1.14) hW' = A(x, h)W := ( Al ( ( ) X) W, W G C 2 , p > 1, 

0 uniformly bounded, with globally separated eigenvalues Ai(.x) = x + i, X 2 = —(x + i). 

Lemma 1.7. There exists T(x,h) on [—L,L] C M, 0 < h < ho, T , T~ l uniformly bounded in C 1 , 
for which W = TZ converts (|1.14l) to a diagonal system hZ' = D(x,h)Z, if and only if 

(1.15) ff x e~ y2 / h ~ 2iy / h 9{y)dy < he~ x2 l h for all \x\ < L. 

Proof. See Section ITT1 □ 


Through condition (I1.15p . we obtain the following counterexamples. 

Corollary 1.8 (Failure of global conjugators). For 0^0 analytic on [—L,L\ x [—i,i\, (11.141) 
possesses a uniformly bounded C l conjugator on [— L, L] as h —» 0 + */ L < 1 and only if L < 1; for 
L = 1, it possesses a uniformly bounded conjugator on [—L,L\ if and only if 6(—i) = 0. 

Corollary 1.9 (Failure of local conjugators for C°° coefficients). Let 9 € C°° be given by 9{x) = 
e~ x forx > 0 and 0 for x < 0, 9 > 0. Then, (11.141) possesses no uniformly bounded C 1 conjugator 
on any interval [—L,L], L > 0. 

Corollaries 11.81 and 11.91 follow in turn from the following estimates proved in Section [3j 


Lemma 1.10. For a ^ 0 analytic on \—L,L\ x [— i,i\, and h —>• 0 + , 

_x 2 

x < he h , 0 < X < L < 1 , 

(116) / e- y2/h ~ 2iy/h a{y)dy ~ h l l 2 a(-i)e~ £ + 0{he~ x2 / h ), \x\ < L = 1, 

/j, 0 + 1 )/ 2 e _ h , 1 < Co < x < L, 

where j is the order of the first nonvanishing derivative of a at z = i. 

Lemma 1.11. For 0 < cq < x < 00 , 9 = e (0, +00), a(y) := e~ y 6 for y > 0 and 0 for y < 0, 

r X 2 /u o- lu 11/0 ~c(s)+d(s)h 1 - 1 / s + 0 (h 2 ( 1 - 1 / s '>) 

(1.17) / e y I iy ^ h a(y)dy ~ h 1 ! 2s e hF~ s , 1 < s < 00, 

J —X 

as h —» 0+, where c(s ) > 0 and Htd(s) = — cos(7r(l — 1/s)) is < 0 for s < 2. 

The symbol a(y) := e~ y ° for y > 0, a{y) = 0 for y < 0 is of Gevrey class & S,T [Lej, |KVl Rmk. 
1.3, p. 3], defined by boundedness of the Gevrey norm ||a|| Sj r := sup^ \dia\(jl) s /T y for some T, 
with s = 1 corresponding to analyticity on a strip of width T about the real axis R. The contrast 
between (11.101) and (11.171) reflects a difference in stationary phase-type estimates for analytic vs. 
C°° symbols a, as quantified in the following more general observation, of interest in its own right. 
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Proposition 1.12. For a E G S,T ° on [-L, L], To, T > 1, |x| < L, and some c = c(Ti,T, s ) > 0, 


(1.18) 


f X e~ y2 l h ~ 2iy / h a(y)dy < /i 1/2 ||a|| T , a e“ c / hl/ ' 
■J —x 


Proposition 1 1.1 2l interpolates between the algebraic 0(h r ) van der Korput bounds for C r symbols 
(roughly, s = oo) and the exponential 0(/i 1//2 e~ 1//ft ') bounds for analytic symbols a obtained by the 
saddlepoint method/analytic stationary phase, as described in Appendix[Bl so far as we know, this 
observation also is new. The lower bounds of Lemma (11.1711 show that (11.1811 is sharp. Specifically, 
for s > 2, = h l ~ 2 / s < 1, yielding ff x e~ y2 l h ~ 2iy l h a(\f)dy ~ hf~ l / 2s e -c(s)/ft. 1/s . for 1 < s < 2, 

the bound is sharp up to the slower-decaying exponential factor 1/s +0(/i 2(1 1/s) ))//j. 1/,s _ 


Remark 1.13 (Failure for analytic projectors). Replacing the diagonal entries of (|1.14ll by Ai = x, 
X 2 = —x, we find by the same argument as for Lemma 11.71 that diagonalization is possible on 
[-L, L] if and only if ff x e~ y2 / h 6(y)dy < he~ x2 ! h for all |x| < L, which clearly fails. Thus, even in 
the case that analytic projectors persist, failure of spectral separation can lead to nonexistence of 
an exact diagonalizing transformation. Note that the radius |x| = 1 of existence of diagonalizing 
transformations in Corollary 11.81 corresponds in the complex plane to the radius at which there 
appears a point z = — i at which Ai = A 2 and separation fails, which is simultaneously a stationary 
point for the phase f = /(A 1 — A 2 ) appearing in (11.1511 . Interestingly, the condition 9(—i) = 0 
determining extensibility up to radius 1 is the condition that A(—i,h ) in (11.141) be diagonalizable. 

Remark 1.14. It is an interesting question whether there holds a general lower bound in (11.181) . in 
which case () 1.18 II would represent an alternative characterization of Gevrey class in terms of the 
F.B.I. transform, analogous to characterizations in terms of the Fourier transform as, e.g., in [FT]. 


1.3. Conclusions. The above results determine what is possible in various settings in terms of 
approximate or exact block-diagonalization, in particular settling the two open problems posed by 
Wasow in his 1965 text |W] on asymptotic ODE whether global exact diagonalization is possible and 
whether analyticity of coefficients is necessary. Namely, we see that, in general, even in the simplest 
case of a system hZ' = A(x)Z + hB(x)Z for which A possesses globally separated eigenvalues, there 
exists (only) a Unite collection of exact locally diagonalizing transformations on neighborhoods 
covering the domain [a, b] of definition of A, B , and this holds in general only for A and B analytic. 

Acknowledgment. Thanks to Gilles Lebeau and Jean-Marc Delort for stimulating conversa¬ 
tions. Thanks to University of Indiana, Bloomington, University of North Carolina Chapel Hill, 
Universities of Paris 7 and 13, ENS Ulm, and the Fondation Sciences Mathematiques de Paris for 
their hospitality during visits in which this research was partially carried out. 


2. Repeated diagonalization and exact local block-diagonalization 

In this section, we compare two methods for obtaining a block-diagonal system from a given 
semiclassical system of ODE, the first approximate, and the second exact. The first can be used to 
compute the second to arbitrary order, while the second gives rigorous validation to the first. The 
approximate block-diagonaization is done globally; the exact block-diagonalization is local. Novel 
aspects of our analysis are the treatment of the point at infinity and of regular-singular points. 

2.1. The method of repeated diagonalization. We start by recalling the method of repeated 
diagonalization as implemented in |MaZ| . by which one may obtain from an approximately block- 
diagonal system with spectrally separated blocks a series of approximately block-diagonal systems of 
successively higher accuracy. For related methods, see, for example, [Ll iFl IW1 [El IBEEKj . Consider 
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an approximately block-diagonal ODE in the semiclassical limit h —>• 0 + : 


( 2 . 1 ) 


hW' = 


M 


o 


0 Al 


li 


22 


Wj + hP 


o e{ 


ol 


Wj, x € [a, b] C 


Proposition 2.1 (Adapted from [MaZj l. Let AP u {x, h), 9j (x, h ) € C r {x ) and C°(h) on [a, 6] x [0, /io], 
uniformly in both coordinates, with r > j > 1, l = 1,2. Suppose, moreover, that A J U have no 
eigenvalues in common. Then, for j < k < r + j, l = 1, 2, and h > 0 sufficiently small, there exists 

Trl h k r k \ 

Me (jr+j-k, conver ti n g (12.11) to 


series of transformations W k = T k W k+ i, T k := 


h k 4 


( 2 . 2 ) 


/i^ = 




Ak 

^22 




Id 


h k e$ 


h k e\ 


w. 


If a = — oo or b = +oo and are exponentially decaying in up to r derivatives as x —>■ Too, 
uniformly in h, then 9 k are exponentially decaying in up to r + j — k derivatives as x —>• Too, 
uniformly in h. If A is analytic on a strip around the real axis and wedges around plus and minus 
infinity, with exponential convergence as Jftx —> Too, then T k , A k , 9 k have these properties as well. 

Proof. We proceed by induction from k = r + 1 up to K , at each step defining T(. such that 

(2.3) 

and setting 

(2.4) 


D k := T k l A k _\T k is block-diagonal, 

A k = T~ x A k _ x T k - h blockdiag{ T' 1 d x T k }, 

Ok = h}~ k ( - T~ l d x T k + blockdiag {T~ l d x T k }, 


has a unique solution of form T k = 


, depending in C 1 fashion on 9. A straighforward 


from which we obtain, evidently, 

To complete the proof, it remains to show that, for 9 bounded and h > 0 sufficiently small, (E3D 

fid 

M Id, 

calculation equating first diagonal, then off-diagonal blocks in the equation T k D k _\ = A k _{T k , 
yields the equivalent system 

a k—1 ^ k ^ k a k—1 | r\k —1 i2 k^k r\k—l„k ukr\k —1 __ k n 

A n a 12 - a 12 A 22 + 0 12 — h a 12 S 21 a 21 - h 0 n a 21 — 0, 

A^ 1 ^ - A k I x + ©M 1 - h 2k a^ - M0M 1 " 


21 


*21"ll 


'21 


*21 v 'l2 


12 


22 


a 


21 


= 0, 


or, written in block vector form in terms of a = ( 012,021 )> 

(2.6) PF{a k , ©fc-iM) := Aa k + Q(a k , 0 fc _i, h) = 0, 

where 


(2.7) 


A 


OL 12 
«21 


- "11 


A^', 1 ai 2 — CI12A: 


A 22 1 ci 2 i — C121A 


fc-i 

22 

k-1 
11 


(2.8) Q(a, 0, /i) = O(|0|)(l + h k \a\ + h 2k \a\ 2 ), Q a {a, 0, h) = O(|0|)(/i fc + h 2k \a\). 

From ra-ra, we have ^7-"(0,0, h) = 0, while, assuming uniform boundedness and uniform 
separation of the spectra of A k f x and A 22 M the decoupled linear operator d a PF(Q,Q,h) = A is 
uniformly invertible, whence we obtain by the Implicit Function Theorem existence of a unique 
small solution a k = Q(Q k -±,h), smooth in both variables. Noting that the properties of uniform 
boundedness and uniform separation of the spectra of A\ x = A k ^ l + 0(h k ) and A 22 = A^f 1 + 0{h k ) 
persist (by smallness of h) throughout the iteration, we are done. □ 




Remark 2.2. The repeated diagonalization expansion may be seen to be a block-diagonal version 
of the classical WKB expansion (11.81) . with p h ’ m analogous to the concatenation 

(2.9) 7m,h '■= T m ■ T m _ i • ... Tj + \ = T m -i,h + 0(h m ). 

By Remark 12.31 the process can be repeated until the system approximately decouples into distinct 
blocks whose eigenvalues all collide: that is, which are either scalar or else possess nontrivial turning 
points. In the case that the initial coefficient matrix has everywhere distinct eigenvalues, the result 
is a complete decomposition into scalar modes analogous to the WKB expansion (11.81) . 

Remark 2.3. A general semiclassical limit problem hW' = A(x,h)W may, by standard spectral 
perturbation theory E, be converted to form (12.11) . j = 1, by an initial block-diagonalizing trans¬ 
formation, W = TW\ such that T~ 1 AT = ^ ^ , yielding (12.11) with A 1 - - = Ajj—h{T~ lr f x )jj , 

e\ = e 2 = 21 , so long as there exist two groups of eigenvalues of A that re¬ 

main separated uniformly in x,h. Specifically, we may take T = {T\,T 2 ), where the columns of Tj 
are bases of the associated total eigenspaces of these two groups, defined by Kato’s ODE 

(2.10) ( d/dx)Tj = [II j,(d/dx)Hj]Tj, 

where Ilj denote the associated eigenprojections and [M,N] := MN — NM the usual matrix 
commutator 0- Thus, there is no loss of generality in starting with the form (12.11) . The projectors 
Ilj, hence the solutions Tj of the linear ODE (12.101) . inherit the same regularity in (x, h) possessed 
by the original coefficient A(x,h). Moreover, if A(-,h) exponentially approaches limits as x —>• 
Too, then Ilj and thus T, converge at the same rate, and 0j < ( d/dx)Tj decay exponentially. 
The transformation T is typically not explicitly computable, but may (similarly as in the WKB 
expansion of Section fl.l.ll) be expressed as a matrix perturbation series in h. 


2.2. Exact analytic local block-diagonalization. Exact diagonalization at a finite point follows 
similarly as in [Wj . as we now describe. 

Proof of Theorem 11.31 Take without loss of generality x* = 0, so that we seek a block-diagonalization 
near x = 0. A straighforward calculation equating first diagonal, then off-diagonal blocks in the 

equation IhT' + TD)Z = ATZ, D = ( 11 +J 1 ^ 11 . ^ „ „ ], yields Ricatti equations 

y U A 22 + nP2'2 J 

holyj = Mia-12 — CK12A22 + @12 ~~ h 2 p ai2©2io:2i — h p Q\\a2i, 

(2.11) ' ,, 

ha 21 = A 22 O 21 — 021^11 + ©21 — /i^a 2 i@i2ai2 — h p Q 22 O 2 I) 


with 

( 2 . 12 ) 


Pi 1 = @11 + /l P 012«21, 

p22 = ©22 + h P Q 21 q 12 - 


Viewed as a block vector equation in a = ( 012 , 021)1 (12.111) has form ha' = Aa + Q(a, ©, h), or 

(2.13) ha' = *4(0)a + ( A(z ) — *4(0))a + Q(a, ©, h), 

where 

a ( a 12\ _ ( ^4llOl2 — 012^22^ 

1 j W ) ~ 1^22021 - aziAuJ ’ 


(2.15) Q(a, 0, h) = O(|0|)(l + h p \a\ + h 2p \a\ 2 ), Q a (a , 0, h) = 0(\@\)(h p + h 2p \a\). 

The eigenvectors of *4(0) may be expressed as tensor products a = \^ l ^ 2 \ and ( 7 ) of 

V 0 j yfoPJ 

eigenvectors of An and A^, an d ^ ie corresponding eigenvalues as the differences in the eigenvalues 








associated with fa, fa and fa, fa. Thus, by separation of eigenvalues of An and A 22 , we find that 
that -4.(0) has no zero eigenvalues. It follows that there is 7 € C, |y| = 1, with argument arbitrarily 
close to zero, for which 7.4(0) has no center subspace. Denoting by II[//ns the unstable/stable 
projectors of 7.4(0), we thus have, for 2 with direction sufficiently close to that of ±7 and rj, C > 0: 


| e h < Ce~ h lri ' ;Rz for 3 ftz > 0, 

e h ~ 1 A(0) z Yi u \ < Ce^ 1 ^, for ftz < 0. 


Defining now II ua = 0 at z* := —M 7 and Il^a = 0 at z* := M 7 for M > 0 real, we obtain by 
Duhamel’s principle the integral fixed-point equation (suppressing dependence of a on h): 


fax) = Tfax) := h 


(2.17) 


+ h 


/ e h - lA ^-^U u ((A(y)-A(0))fay) + Q(a,Q,h)(y)) dy 

J 2 (* 

1 f e h - lA ^ x -^U s {(A(y) - A(0)) fay) + Qfa, ©, h)(y)) dy, 

J z* 


defined on the diamond 

(2.18) V :={x : | arg ((x - z*)/ 7 )!, I arg ((z* - x)/fa \ < e} 

containing a neighborhood of the origin, with 0 < e,M <C 1. Noting that [x — y) has angle 
arbitrarily close to that of ±7 for x E T> and y G \z*,x\, [x,z*\ for e > 0 sufficiently small, we have 
by MD the bounds | e fc _ M(°)(*-i/)n s | < Ce~ h ~A^-y) and | e ft-M(o)(*-f/) nt/ | < c^-^x-y) in 
(I2.17p . whence we obtain by integrability of h~ 1 e~ h r]t over t G M + the estimates 
||T(a)|| L cx>(x>) < C\\A — ^l(0)||x,oo( :D )||a|| i oo(X)) 

+ C'II©I|l«(d) + h 2p \\ a \\ 2 L°°(v)^ > 

||T(ai — ot2)\\L°°(v) < C , (||“4 — A (0)|| loo (xi) + ||©||l°°(x>) (h p + h~ p sup || “ill L°°(X>))) 


(2.19) 


X ||«1 — a2111,00 ( 23 ). 

Taking M, hence |A — -4(0)|, and h sufficiently small, we find from (|2. 191) 41 that T takes the ball 
S(0,2C||©|| L °°(v)) t° itself, and from (I2.19D (ii) that T is contractive on R(0, 2Cj|0||£°°(£>)) with 
contraction constant < 1/2. It follows that (12.171) determines a unique solution for M,e <C 1, 
which, moreover, is bounded as claimed. Regularity with respect to parameters is inherited as 
usual through the fixed-point construction via the Implicit Function Theorem. □ 


2.3. Exact block-diagonalization at infinity. Exact diagonalization at infinity follows similarly, 
but with an important modification having to do with restriction to a wedge. The argument fails 
for functions merely analytic on a strip; analyticity on a full wedge is needed in the proof. 


Proof of Theorem\1.4\ Under assumptions (11.120 . consider the equation 


( 2 . 20 ) 


, dW 

h ib = 


An 

0 


0 


A 


22 


(x,h)W + h p O(x,h)W. 


Following the same steps as in the proof of Theorem 11.31 we obtain as in (1 2.13 [) equation 
(2.21) ha' = . 4 ( 00 ) 0 : + 0(\A(z) — .4(oo|)o + Q(a, 0) 

h p an 


for entries a = ( 072 , 0 : 21 ) of a diagonalizing transformation T = 


, A as in (12.141) . 


/ 

Ji p a2i I 

Again, we note that the eigenvalues of .4(oo), formed by differences in eigenvalues of An ( 00 ) and 
^ 22 ( 00 ), do not include the value zero. Choosing e > 0 so that there are no eigenvalues of .4(oo) on 
the rays with angle ±(vr/2 + e), we can thus divide the eigenvalues of 7l(oo) among three subsets 
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in a way that that persists under small variations in h (or the suppressed parameter q of ( 11 . 11 ) ). 
namely: I. eigenvalues with argument lying strictly between ir/2 + e and 3n/2 — e; II. eigenvalues 
with argument strictly between —e and 7 r /2 + e; and III. eigenvalues with argument strictly between 
— 7 r /2 — e and s, with associated projectors II/, II//, and II///. 

Taking now e > 0 sufficiently small, and defining the wedge 

( 2 . 22 ) Wm', e : {x : — e < arg(x — M') < e}, 

and direction vectors j± := e ±l£ , and taking M' > 0 sufficiently large, we can, similarly as in the 
finite case, obtain a as the unique solution in L°°(Wm', e ), of the fixed-point equation 

(2.23) 

a(x) = T'a(x) := h~ l f e h — A(oo)a)) + Q(a,@,h)(y)) dy 

JM' 

+ h - 1 f e /.-U(oo)(x- y )n n ((A(y) - A(oo)a) + Q(a , 0, h){y)) dy 

J x+ 7 _ (+oo) 

+ h~ l f X e h - lA ^ x -^U in ((A(y) - ^l(oo)a) + Q(a, 0, h){y)) dy , 

J £+7 + (+oo) 

where the contour integrals are understood to be along straight lines. 

For, noting that multiplication of A(oo) by y± rotates its spectrum by angle ±£, and recalling 

that eigenvalue in region II have angle strictly between — e and 7t/2 + e, we have that 7 _4(oo)II// 

restricted to range of II// has eigenvalues of angle strictly between —2e and 7 r /2 and, similarly, 
7 + *4(oo)II/// restricted to range of II/// has eigenvalues of angle strictly between — 7t/2 and 2s, 
whence for e less than n/2 both have eigenvalues of strictly postive real part. Likewise, for — s < 
9 < e, e‘ e A(co)Ai restricted to range of 11/ has eigenvalues of strictly negative real part. It follows 
that, analogously to (12.161) . 

| e h-V e M( 0 ) n /| < Ce -fc-V f or | 0 | < e , t > 0 real, 

(2.24) \e h ~ lry + tA ^U n \ < Ce h ~ lyt , for t < 0 real, 

\e h ~^- tA{ - 0) U in \ < Ce h ~ lljt , for t < 0 real, 

for some C,r) > 0, yielding analogously to ()2.19l) 

(2.25) 

||r(a)|| L oc( WM , e ) < C\\A - ^4(00)||£,oo(w M/ie )||Q:||Loo(yv M , ]e ) 

+ C'||0|| L oo (V v M , |e )(l + h p \\a\\ L °°(y VM ,^ + h 2p \\a\\ 2 L0O ( WMi ^, 

\\T(ai - a 2 )||//x>(w M , >E ) < c(\\A - ^(oo)|| L oo( W m; j + ||0||loo(w m/jE ) (h p + h 2p sup ||aj||/,oc.(w M;ie ))) 

X ||«i - CK2IIL°°(VV M , e )- 

Taking h sufficiently small and M' sufficiently large, hence \A — ^l(oo)| sufficiently small, we hnd 
from (12.250 (41 that T takes the ball B{Q,2C\\Q\\ L oo^ m , to itself, and from ()2.25jl (ii) that T 
is contractive on S(0, 2C'||0||/,c»(yv M )) with contraction constant < 1/2. It follows that (12.170 
determines a unique solution for M, £< 1, which, moreover, is bounded as claimed. Regularity with 
respect to parameters is inherited through the fixed-point construction. Finally, straightforward 
estimates using exponential convergence (decay) of 0 in y, and exponential decay of | A(y) — »4(oo)| 
and propagators e h in |3ft(x — y)\ yield exponential convergence (decay) of a. □ 

Remark 2.4. The introduction of oblique contours [x, x+ 7 ±(+oo)], made possible by the assumption 
of analyticity on a wedge, is what makes possible the subdivision of eigenvalues of *4.(oo) yielding 
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strict exponential decay. With analyticity only on a strip, we would have a problem with algebraic 
growth in case that *4.(oo) had a pure imaginary eigenvalue with nontrivial Jordan block. Less 
critically, we would have at best neutral decay of propagators, hence would lose a power of h in the 
exact conjugation, since h~ 1 e~ r,t / h is integrable for t € M + only for 3 £77 > 0 (see [Zll ILWZIj L 

Remark 2.5. By 0(e~ ri ^ x ~ y ^) decay of propagators, the contribution to (|2.23p from \x—y\ > \x\/C 
is negligible for any C > 0, whence, taking C > 1, we recover from decay (convergence) of 6 decay 
(convergence) of a, whether or not the rate of decay (convergence) is exponential. In particular, 
for algebraic decay (convergence) of 8, we obtain decay (convergence) of a at the same rate. 

2.4. C r version in the case of a spectral gap. By a modification of the “variable-coefficient 
conjugation lemma,” Lemma 4.3 of jZl] (see also Remark 7.2 of the same reference), we may 
alternatively obtain existence of a block-diagonalizing conjugator for Ajj, 0 real-valued and merely 
C r , r > 1, again on [M, 00 ), M » 1, under the additional assumptions that (i) 3fL4n > JM 22 or 
vice versa, where IRN := (1/2 )(N + N *) refers to the symmetric part of a matrix N, and (ii) 0 is 
exponentially decaying as x —>• + 00 , uniformly in h. If p > 2, we may, further, take M with any 
value (even egative), for h > 0 sufficiently small, where, Ajj and p are as in (II.111) . (Recall, Remark 
El that exponential decay condition (ii) may be arranged so long as there is spectral separation 
(as opposed to gap) between blocks.) The proof amounts to the observation that we may in this 
case use the fixed-point construction (12.171) restricted to x € R and z* := — M, z* = M, noting 
that (12.161) hold with 7 = 1 and rj = 0. This clarifies perhaps what can already be done with C r 
coefficients and what is gained from analyticity. We show in Section 0] that these distinctions are 
actual and not only apparent; that is, the hypotheses for analytic vs. C r versions are sharp. 

Remark 2.6. If an analytic function has an accumulation point of purely imaginary values on the 
real line then it is imaginary-valued on the real line (inspection of coefficients at the accumulation 
point). Thus, crossing of real parts are isolated, and one can use the real-valued variable coefficient 
conjugation lemma to treat intervals between, so long as the neutral numerical range condition 
JiAii > JM 22 is maintained- roughly, so long as neutral eigenvalues stay semisimple. An advantage 
of this approach is to know the size of intervals on which conjugators exist. On the other hand, the 
C r conjugators so constructed are not in general analytic, even if coefficients Ajj, 0 are analytic, 
unless the neutral numerical range condition extends to a complex strip. 

2.5. Diagonalization at a (finite) singular point. 

Proof of Corollary \1.5l This case may be converted to the problem treated previously in Theorem 
[L4l In coordinates x := — In 2 , we obtain 

(2.26) h ^r={ A 0 + 

Aj, © analytic in x on a half-plane > M, and convergent as Jix —» +00 (z —>• 0). Applying 
Theorem 11.41 we obtain a conjugator that is analytic in x = — In 2 ; on a wedge x € Wm',£ as in 
(I2.22p for M 1 > 0 sufficiently large and e > 0 sufficiently small and convergent as —>• 00, hence 
analytic in z on a slit disk about z = 0 with cut along the negative real axis and continuous at 
z = 0. □ 


3. Oscillatory integral estimates 

As explored further in Section f4.ll existence of conjugators is related to decay rates for certain 
oscillatory integrals. For later use, we carry out here some estimates needed for our analysis. 
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3.1. Variation on a classical sum. We first generalize the estimate j^ 00 e \ e h V dy r^j h}/ 2 e- i/ fc , 
obtained by direct Fourier transform to general analytic symbol a(y) and bounded domains. 

Proof of Lemma \l.lQ[ Consider ff e _ We — a*a(y)dy, where a is holomorphic on [—L,L] x [— i,i\. 
For x > Co > 1, we obtain the estimate by the Analytic Stationary Phase Lemma IB . 1 1 together with 
Remark lB.2l ih defining the contour C = [x, — i T e] U [—i T e, — i — e] U [—i — e, —x], 0<£<1 real, 
and noting that [—x, x] U C forms a closed contour in the domain of analyticity [— L, L\ x [— i, i] of 
a, on which the phase (j>(y) := —y 2 — 2 iy has a single quadratically nondegenerate critical point at 
y = —i and, for e sufficiently small, satisfies %t<j)(y) < 9i0(— i) = — with 3i(/>(a), 3?0(6) < (R</>(—*). 

The estimate for x < L < 1 follows, similarly, defining the contour C' = [x, — i|x| + e] U [—i|x| + 
£, —z|x| — e] U [— i\x\ — £,—x], and observing that the phase (f> has no critical points on C' with 
5R <j)(y ) < e~ x ! h throughout. Expressing = f c , using Cauchy’s Theorem, we obtain by a 
nonstationary phase computation, integrating by parts, finally, 

(e-^/ h y a { Z )/(/)'( Z ) dz 

= he-^ z )l h a( Z )/<t>'( Z )\ x _ x + h [ e-^ 2 )/ /^ (a(2)/^(z)) , d2</^e- a;2 / ,l . 

Jc 

For x < 1 = L, using the original contour C, a nonstationary phase estimate like (13.11) shows that 
the contribution from the parts [x,—i + e] and [—i — £,—x] is < he~ x ! h , while the contribution 
from [—i + £, — i — £] may be estimated directly from (IB. II) as ~ (L 1 / 2 a(— i) + 0{h))e 1 / h , whence, 
summing, we obtain ff x e~^ z ^ h a( Z )d Z r^j /i 1 / 2 a(— i)e 1 / h + 0{he x ! h ) as claimed. □ 


[ e~^ z)/h a( Z )d Z = -h [ 

Jc Jc 


3.2. A C°° oscillatory integral. We next investigate the integral 


[ e V h 2l ha{y)dy= [ 
J-x Jo 


= e h 


la (y)dy 


with symbol a(y) = 0 for y < 0, a(y) = e y for y > 0 that is C°° but not analytic. 


Proof of Lemma U.lR Take without loss of generality x = +oo, noting that the resulting error is 

V 2 

negligible, denoting the resulting quantity 1(h) := f 0 °° e h 2l ha(y)dy. Define 

i7r(l —1/ s) 

(3.2) a = 1 - 1/s € (0,1), j3 = e -a-. 


Case (i) (1 < s < 2). Deforming the contour [0, +oo] using Cauchy’s Theorem to the contour 

z = h a f3t, t £ (0, Too), 

we obtain 


(3.3) 


1(h) ~h a e 

Jo 


— — t ^ J r i[3h a t^) 


h l-° 


'dt. 


where the real part of the phase, ?R,i(3(—2t — t~ e T if3h a t 2 ), has a unique, quadratically degenerate 
maximum for h = 0 at to = 2 ~C- 1 / S ). Applying the Analytic Complex Stationary Phase Lemma, 
Lemma IB. 11 more precisely the generalization of Remark IB. 21 3) to an h-dependent phasejl we 


obtain finally 1(h) ~ h a h 2 e 5R b/ 5 )( 


-2t 0 -t~ a +i^h ol +0(h^ a ) 


h 1 - 01 ), yielding the result. 


Case (ii) (s > 2). In this case, 3?(— /3 2 ) = — cos(7r(l — 1/s)) > 0, and so we cannot move the 
contour as in the previous case to a ray in direction j3\ indeed, the integral in (13.31) is not convergent. 
Instead, we first move to the furthest possible ray z = h a e~ m4 t, t £ (0, Too) on which the infinite 


2 With also C = [b, a], in the notation of the lemma. 
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integral converges (and for which the contribution at infinity vanishes for all rays between, justifying 


the shift of contour), obtaining 1(h) h«Ce 


ftl-o 


" dt for /3* = e m4 , then truncate to 


(3.4) 


1(h) ~ K 


f 


ip*(-2t-t~ u +ip J ,h°‘f s ) 
-> J.1 — a 


'dt 


for L 1 at the cost of a negligible 0(e cL ! h s ) contribution. Finally, we estimate (13.41) using the 
/i-dependent stationary phase lemma, with countour C = \/3*L,(3L\ U \/3L, 0], /3 as in (13.21) . yielding 

' -2t 0 -t~ 9 + iph a +0(h 2a ) ^ 


h 1 - 0 


as in the previous case. 


□ 


the same estimate 1(h) ~ h a h 2 ^ l h)(- 

Remark 3.1. The rate of exponential decay obtained in Lemma [LIT] may be recognized as e^( z * LFL. 
where z*(h) is a critical point of the augmented phase 4'(y, h) := (—y 2 /h — 2 iy/h — 1 /y e ) obtained 
by including the symbol e~ l ! y as part of the phase. 

3.3. Gevrey-regularity complex stationary phase. 

Proof of Proposition \1.12\ First observe that, by estimate (11.171) . we may construct cutoff functions 
£ of Gevrey (s, T) class for arbitrary s and (by rescaling) T, hence, muliplying a by such a £ 
supported near the origin and periodically extending, the observing that the error incurred is 
negligible, we may restrict to the case of a periodic. For periodic functions Foias and Temam have 

derived an equivalent version ||a||* T := \Jj2jez( 1 + |j|) 2 e 2:z ^ 1/s \dj\ 2 of the Gevrey (s,T)-norm 

[FT], [EB (2.4), p. 4], where dj denotes Fourier transform, with T in the analytic case s = 1 
corresponding to width of the strip of analyticity about the real axis R. For T > 1, and ||a,||* 

finite, the complex stationary phase integral L a (h) := f R e~ y2 / h ~ 2w / h a(y)dy may be seen to satisfy 
L a (h) < C||a||^ s h ka e~ c / hl/a , c, C > 0, by triangulation with the function a, analytic on a strip 
of width T > 1, obtained by Fourier truncation of a at modes |£| < 1/Th , using the analytic 
stationary phase estimate 

L a (h) < < /i 1 / 2 |a|| 7 . 0 iS e~ 1 ' /,,1/s 

following from (11.161) . the fact that sup|a| on the strip [—L,L\ x [— i,i\ of width 1 is bounded 
by Hall^-p and the straightforward estimate ||a||y 1 < IMI 7 1 s e 1 / /l ~ T ’°( T /L 1/s , together with the 
truncation error bound 


sup 

cce[—L,L]cK 


— a I < 


a — a 


E 

\j\>T/h 


N < 


/ E a + lil) 

\j\>T/h 



E a + ui )- 2 

\0\>T/h 


< ll«llI.ryFPFe-™ 1 '- < ft'/ 2 ||a||; iT e-V». 2 '..)/'.'", 

where we have used in the first inequality Hausdorff-Young’s inequality, in the second Cauchy- 
Schwarz’ inequality, and in the third the definition of || • ||* T . The latter bound gives evidently 

\L a (h) - L a (h )| < sup \a — a\ < /i 1/ -||a||* T e _c(To ’ T ’ s)/hl/s , 

aiGM 

from which we obtain finally that, for some c > 0, 


\La(h)\ < \L a (h)\ + I L a (h) - L a (h )I < h l / 2 \\a\\l To e- c ' hl/ 


□ 


Remark 3.2. The Fourier truncation argument above gives a general way of interpolating between 
results for analytic and C r functions. For example applied to the analytic interpolation bound 
interpolation error < C'||/||i i t(T+v / 1 + T 2 )~ N of |DYj (adapted from [T]). with cutoff |^| > N, this 
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yields the more general bound error < C||/|| ^(T + y/l + T 2 ) Nl/s for Chebyshev interpolation 
of periodic functions / on [—1,1] with N mesh points. 

Remark 3.3 ( C r stationary phase). Similarly, the W r, °° symbol a(y) = 0 for y < 0, a(y) = y r for 
y > 0 yields by a standard nonstationary-phase argument involving r integrations by parts 

/ X PX 

e ~ y2 l h ~ 2iy / h a{y)dy = / e~ y2 ^ h ~ 2iy ^ h a(y)dy ~ h r for 0 < Co < x, 

-x Jo 

verifying sharpness of the nonstationary van der Korput bound for W r)CO symbols a [M] , 

4. Counterexamples 

We complete our analysis by carrying out the counterexamples described in the introduction. 

4.1. Block-diagonalization and decay of oscillatory integrals. Recall the triangular system 


hW' = AW := 


M 

0 


h p e 

A2 


w 


introduced in (|1.14D of the introduction, with Ai(x) = x + i, A 2 = — (x + i), W € C 2 , p > 1. 

Proof of Lemma We seek a coordinate change W = TZ with T, T 1 uniformly bounded in 
C 1 , converting (11.141) to an exactly diagonal system HZ' = DZ. From p > 1, uniform boundedness 

in C 1 , and the relation D = T~ 1 AT — hT~^T', we find, first, that D = (^ . 

\ 0 A 2 + 0(h)^ 

Comparing 0(1) terms in the off-diagonal entries of T~ 1 AT } we have also that Tn7i2> T 21 T 22 ^ h, 
and, since detT must be bounded below by boundedness of T” 1 , T 11 T 22 — T 12 T 21 > 1. 

The latter observations imply that either (i) Tn, T 22 ~ 1 and T 12 , T 21 h, or (ii) T 12 , T 21 ~ 1 and 

Tj i, T ‘22 JS h. Noting that case (ii) may be reduced to case (i) by exchanging T for T = T J 

we may assume without loss of generality (i). By a further, diagonal transformation rescaling the 

diagonal entries of T to one, we may arrange that T = , with (a, a)(x, h) < 1. Comparing 

off-diagonal coefficients of T~ 1 AT — hT~ 1 T ', we find by an induction argument, finally, that 


(4.1) 


T = 


1 hPa 
hPa 1 


D = 


Ai + 0(h p ) 0 

0 A 2 + O(hP) 


(a,d)(x,h) < 1. 


(Alternatively, for purposes of this argument, we could take without loss of generality p = 1.) 

This reduces us to the case treated in Theorem ll.3l with @n = ©22 = ©21 = 0, ©12 = 0, a \2 = a 
and Q! 2 i = a, in which (12.111) specializes to 

ha = (Ai — A 2 )a + 9 , 

(4.2) 

ha = (A 2 — Ai)d — h p a9a, 

with d\ = Ai + h p 9a and ^2 = ^ 2 - Noting that ()4.2|1 (i) is independent of (14. 2D (ii). and that (|4.2jl (ii) 
is consistent with a = 0, we find that there exists a diagonalizing transformation if and only if there 
exists a triangular one, a = 0, which exists if and only if there is a solution a of (12.lip (i) that is 
uniformly bounded as h —>• 0 -1 ", equivalent by Duhamel’s principle to uniform boundedness of 


(4.3) 


a(x, h) = h~ l [ e h 1 ^ (Xl ~ X2){z)dz 9(y)dy - e h 1 fo(^-^)dz a ^ ^ 

Jo 

= e ( x2 + 2ix )/ h (h~ 1 J* e -( y2+2iy V h 6(y)dy - a(0, hfj . 
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(<=): Noting that e 2lx a(x ) — e 2lx a{—x) = e x2 / h f* x e ^ y2+2w ^ h 6(y)dy, we see that (14.31) is 
uniformly bounded as h —» 0 + only if J// x e~^ y2+2iy ^ h 9(y)dy < he x2 ! h . 

(=>): Choosing a(0, h) = h~ l e~^ y2+2iy ^/ h 9(y)dy, we obtain from (14.31) 

,-L 

(4.4) a(x,h)=e( x2+2ix V h h- 1 e -fe 2 + 2 *v)/H(y)dy, 

J X 

or, alternatively, 


(4.5) 


a(x, h) = e (x +**)/ h h~ 1 / e~ (y +2iy) / h 9{y)dy 

J X 

+ e ^ 2 +^)/h h ~l j L e -(2/ 2 +2 iv)/h0(y)dy. 


For 0 < x < L, a nonstationary phase estimate as in the proof of Lemma 1 1.101 (case \x\ < 1) yields 
uniform boundedness of the righthand side of (14.41) as h —» 0 + ; for — L < x < 0, the same estimate 
yields uniform boundedness of the f x term on the righthand side of (14.51) . while assumption (11.151) 
yields uniform boundedness of the term. Thus, a is uniformly bounded for all |a;| < L. 

□ 


Remark 4.1. From the case A j constant, we see that in general o(0, h) 7^ 0, or T(0, h ) 7^ Id. 

4.2. Global counterexample. 

Proof of Corollary ll.ffl For any nonvanishing 6 analytic on [—1,1] x [— i,i\ and continuous on 
[— L,L\, we hnd by Lemma 11.101 that condition 11.151 is violated for L > 1. □ 


4.3. C°° counterexample. 


Proof of Corollary ] 1.81 For the Gevrey class s function defined by 6{y) = 0 for y < 0 and 6{y) = 
e~ y for y > 0, 9 = —k-, we hnd by Lemma 1 1.11 1 that condition II. 151 is violated for any L > 0. □ 

Remark 4.2 ( C r counterexample). For the W r ’°° function defined by 9{y) = 0 for y < 0 and 
6{y) = y r for y > 0, we hnd by Remark 13.31 that condition 11.151 is violated for any L > 0, giving a 
particularly elementary counterexample to existence of diagonalizers in the G r_1 case, r arbitrary. 


4.4. Regular-singular point counterexample. Finally, for singular ODE w' = (1 /z)A(z,h)w, 
with A analytic on the disk R(0,r), we show that there may not always exist an analytic block- 
diagonalizing conjugator on the whole disk, but only on a slit disk as described in Corollary 11.51 


Example 4.3. Consider hzW' = ( q q ) ID + h 


0 zcj) 

0 0 


W, f> analytic. Looking without loss of 


generality (following Section 14.11) for a triangular conjugator T = 


1 ha 
0 1 


, we obtain the ODE 


ha' = a/z + (f. Assuming Taylor expansions a(z;h ) = Ylj a j{h)z :) , cj>(z) = Ylj ■> we find, 
comparing coefficients of like order, that no = 0 and (jh — 1 )aj = <pj -1 for j > 1. Thus, there is a 
formal solution for all 0 < h < ho if and only if 1 f> is polynomial, in which case there is a polynomial 
solution a. For, otherwise, taking hj = 1/j, we have a sequence hj —>• 0 + for which there is no 
solution whenever <fj—i 7 ^ 0. Rewriting the ODE as (z~ha)' = hz~^ tpjz 3 and integrating term 
by term, we see explicitly the appearance of a (fj-iz 3 log z term when h = 1/j and 4>j-i / 0. Thus, 
in general there may not exist an analytic diagonalizer for the regular singular point problem (11.131) . 
On the other hand, if blocks An (0,0) and ^ 22 ( 0 ,0) have no eigenvalues differing by a real value, 
one can show that an analytic solution ( 011 , 0 : 2 ) exists for \z\ < r, 0 < h < ho, by solving for the 
formal power series as above, and showing convergence by direct estimates of coefficients. 
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Appendix A. Analyticity of traveling-wave profiles 
A. l. An Analytic Stable Manifold Theorem. Consider a complex analytic ODE 
(A.l) (d/dt)u = f(u), t € C 1 , u € C n , , / € C u : C n —>■ C n , 

defined in a neighborhood of an equilibrium u*, /(it*) = 0, with associated linearized equation 
(A.2) v' = Au, A := df(u *). 

Associated with A, define the stable subspace as the direct sum of all eigenspaces of A associated 
to stable eigenvalues, or eigenvalues with strictly negative real part. Likewise, define the center, and 
unstable subspaces E c and as the direct sum of eigenspaces associated to neutral and unstable 
eigenvalues, respectively, i.e., eigenvalues with zero and strictly positive real part. This gives a 
decomposition C n = E s © E c © E u of C n into subspaces that are invariant under the flow of (IA.2I) . 

Defining associated (total) eigenprojections II S , II C , and LU as the sum of all eigenprojections 
associated with stable, neutral, and unstable eigenvalues, respectively, we have the standard bounds 

\e At U s \ < C(r ? )e“ r?|m| , Rt > 0, \Qt\ < v\Rt\, 

(A.3) \e At Ii c \ <C(q)e m \%t\ < v\Rt\, 

\e At u u \ < c(9)e~^ m i, Rt < o, \st\ < v\Rt\, 

for any r/ > 0 strictly smaller than the minimum absolute value of the real parts of stable and 
unstable eigenvalues, 6 > 0 arbitrarily small, and some v > 0. 

Theorem A.l (Analytic Stable Manifold Theorem). For f E , there exists local to it* a 
stable manifold M. s tangent at u * to E SJ expressible in w := u —u* as a C u graph : E s —>• E C ©E U; 
that is (locally) invariant under the flow of (I A. ID and uniquely determined by the property that 
solutions in A4 S decay exponentially to u* in “approximately forward” time, in the sense that 
|v;(t)| < C(fj)e~^ t )w(0)\ for Rt > 0, |3#| < n\Rt\, and any 0 < fj < rj, rj, v as in da. 

Proof. Defining w := it —it*, we obtain the nonlinear perturbation equation w' = Aw + N(w), where 
A := (df / du)(ufl) is constant and N(w) := /(it* + w) — /(it*) — {df / du)(u*))w is the first-order 
Taylor remainder for /(it) — /(it*), satisfying N( 0) = 0, ( dN/dw)(0 ) = 0, hence 

(A.4) lim | N(w)\/\w\ = 0, lim | (dN/dw)(w)\ = 0. 

|«j|-> 0+ |iu|—> 0 + 

Applying projections Ilj, j = s,u,c, and using the fact that these commute with A, we may 
coordinatize as three coupled equations (Iljit;)' = A(Hjw) + HjN(w) in variables w s := II s it>, 
w c := n c it;, and w u := II u w. Considering IljiV as inhomogeneous source terms and applying the 
variation of constants formula, we obtain 

(A.5) Ujw(t) = + l e^ (t “ s) n i A(s, itt(s)) ds, 

J toj 

j = s, c, it, so long as the solution w exists. 

Under assumption \w(t)\ < Ce~^ 1 for Rt > 0, \$st\ < n\Rt\, we find by (1A.3I) that \e A ^ t ~ r 'l'n.jw(T)\ 
decays exponentially as Rt —>• +00 for j = c, u with |9(r — £)| < v\R(t—1)\ and Rt > 0, |9?i| < v\Rt\, 
while e A ^~ s ^IijN(s, w(s)) ds (since |A(ic,s)| < Cjie(s)| for |u;| uniformly bounded) converges 
to a limit. Thus, taking foj —> +00 for j = c,it, in (1A.5|) . we find that the first (linear) term on 
the righthand side disappears, leaving Hjw(t) = — e A ^~ s ^IJjN(s, w(s)) ds. Choosing to, s = 0 
and summing the three equations (IA.5I) . we obtain for w s := n s iu(0) the integral fixed-point 
representation 

rt /*+00 

(A.6) w(t) = T(w s , w) := e At w s + ® e A ^ s ^H s N(w(s)) ds — (h e A ^~ s ^H cu N(w(s)) ds, 
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valid for all solutions decaying exponentially in approximately forward time at rate Ce -r,5Rf , where 
n cu := II C + n u denotes the total eigenprojection of A onto the center-unstable subspace E cu := 
E C ©S U . Here, dt denotes any contour integral lying within the cone $t(s — t) > 0, |9?(s —f)| < 
n|9(s — t )| with Jis —> +oo. 

Define the time-weighted norm ||/||jj := sup Kt > 0j e^ m \f{t)\, noting |/(f)| < e~ vm ||/||^ 

for Jif > 0. \Cst\ < With respect to this norm, we find by a straightforward (and standard, 

in the real setting) computation using (1A.3[) that, for |w s | and <5 > 0 sufficiently small, the integral 
operator T is Lipschitz in w s and contractive in w from the || • ||^-ball B(0,5) to itself: 

(A.7) \\T(w s ,wi) - T(w s ,w 2 )\\rj < {l/2)\wi - w 2 \\ij. 

Thus, by the Contraction Mapping Principle there exists a unique solution w = w(w s ) £ B(0,5) 
in Bff := {/ : \\f\\fj < oo}. Appealing, finally, to the analytic dependence of analytic fixed point 
mappings on parameters, we obtain that w(w s ) is a C u function from E s C R n —>• Bf r 

Evidently (computing n s u>(0) = w s by applying n s to the right-hand side of (1A.7I) and using 
n s n cu = 0 ), the function w(w s ) takes each point in T, s to a unique exponentially decaying solution 
on Jit > 0, \Qt\ < v\iRt\ with initial data satisfying n s w;( 0 ) = w s . Thus, defining < f > (-u; s ) := 
n cu u;(r/; s )|t = o, w e obtain the desired C u map, whose graph, a C u manifold, consists, locally, of 
the set of data with exponentially decaying solutions, or, equivalently, the invariant set of orbits 
decaying exponentially in approximately forward time. Tangency, $(0) = 0 and cid^O) = 0, then 
follow by uniqueness together with the fact that A"(0) = 0, (dN/dw)( 0) =0. □ 

Remark A.2. The above argument can be recognized as the standard fixed-point construction for 
ODE on the real line, together with Cauchy’s Theorem plus the observation that the standard 
linearized bounds (|A.3D hold on a wedge. In the absence of a center subspace, essentially the same 
argument yields that the stable manifold is uniquely decribed as the set of solutions merely bounded 
in approximately forward time. Indeed, setting u = 0, we find that all solutions merely bounded 
and close to w* in forward time for t real must lie in the stable manifold as well. 

A.2. Application: analyticity of profiles. Applying the above stable (unstable) manifold the¬ 
orem to a profile ODE, we obtain immediately a result of analyticity on an approximately forward 
(approximately backward) sector. Consider a general profile ODE 

(A. 8 ) v! = f(u), 

written as a first-order system, with / analytic on the set «£Wc C n in question and ' denoting 
(d/dx ). Here, the connection problem can be either a heteroclinic/homoclinic one on the whole real 
line, or a boundary-value problem on a forward or backward real half-line. In any case, provided the 
desired endstate in forward real time is a hyperbolic rest point for the ODE restricted to the real 
axis, we find from Remark IA.2I that any connecting profile (since bounded and close to n*, at the 
least) must belong to the complex analytic stable manifold described above, consisting entirely of 
orbits analytic and exponentially decaying in an approximately forward sector in x. More generally, 
this holds for any exponentially decaying profile, even in the presence of a center subspace at li*. 

Appendix B. Complex saddlepoint estimate/method of stationary phase 

Lemma B.l (Analytic Stationary Phase Lemma: adapted from [Si IPWj ). Let C U [a, b] C C be an 
oriented continuous closed curve in the complex plane, a, b real, with standard orientation on [a,b\, 
enclosing domain Ll, with cj), A analytic in x on Cl and A Lipschitz continuous in h for 0 < h < ho, 
ho > 0. Suppose further that there is a quadratically nondegenerate critical point zo € C mt of <f> 
such that Ht(j) has a nonstrict maximum on C at zo, with C described by a C 2 curve 7 (f), t € [— e, e] 
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in the vicinity of zq, 7(0) = zq, and this is the only critical point on C. Then, 


(B.l) 1(h) := JI e ^ A(x,h)dx = h}^ 2 e h (A(zq,0)^ 

as h —» 0 + , where the square root is chosen so that — 7 , (0) and 


<K~ 0 ) 


2vr 


-&'(* 0 ) 


+ 0(h~ 1 / 2 )) 


lie in the same half-plane. 




Proof. By the Cauchy integral formula, f CU [ a{) ] e~h~ A(z,h)dz = 0, whence, rearranging, 

7~)-<7z n) 


I(h)e h = — 


L 


e h A(z, h)dz, 


<p( z ) — <t>(z 0 ) _ _ __ 

where, by assumption, e h is uniformly bounded. Substituting A(z, 0) for A(z, h) gives a 

negligible 0(h) error, reducing to the classical case treated in | 51 IPW|. Next, loosely following 
[FW] , deform 7(f) to 7(f) := 7 (t) — ?yx(f)V(!R 0 ), where x(t) is a C°° cutoff function vanishing at 
t = =t e and identically 1 in a neighborhood of zq, and 77 > 0 is taken sufficiently small that the image 
of 7, and the domain between the images of 7 and 7, lies in the region of analyticity of A(-, 0 ), <f: 
that is, an approximate gradient descent. It is readily seen that, for 77 sufficiently small, 7(0) = zq 
and $t(j)( 7) has a quadratically nondegenerate maximum at t = 0, while $t(<fi(y(t)) — (j>(zo) < 0 for 
t € [— e,e\. Now, for \t\ > h 5 / 12 , we may write the integral over 7(f) as 


/ 


<K Z )-<H Z n) , . . 

e h A(z, 0) = 


4>(z) — <t>(zri) 

hd z (e - h a ~)A(z, ti)/4> (z) 


and obtain a uniform 0(h) bound by integrating by parts. For \t\ < /i' 5//12 , meanwhile, we have 
t 3 /h « 1, so that e°d 3 / h ) = 1 + 0(t 3 /h). Taylor expanding, therefore, and keeping only 


n 4>(z)—4 >( z q) _ _ __ _ 

lower order terms, we hnd that on this region — J e h A(z,0)dz reduces modulo O(h) errors 
to the exact integral f R e * ( 


■A(zq, 0)(—7 / (0))dt = h 1 / 2 A(z 0 ,0)^ , giving the result. 

r]t 2 rjt 2 

Here, we are bounding remainder terms using (|f| + \t 3 \/h)dt < /i -1 / 2 f e - "^" dt = 0(h~ 3 ) 


nh 10 / 12 


and fi t i >h s/i 2 e 2h dt = 0(h}^ 2 e h ) <C h x , and using the change of coordinates z — zq = 
re 4uce the exact integral to a standard Gaussian integral 

V'(*o)l f 

J R 


H(^ 0 ,0) 1 


-<t>"(zo) 


\cj>"(0)\w 2 

e zh dw = A(zq, 0L 


2tt h 


(/)"(zq) ' 


□ 


Remarks B.2. 1. Lemma IB.II mav be generalized in straightforward manner to the case that C 
contains multiple critical points, or critical points at its endpoints [ PW] , yielding an asymptotic 
expansion as the sum of contributions of each critical point (with endpoints counting half). Likewise, 
for A(z 0 ) = 0, and making the additional assumption of strict inequality 3?0(a), 3?^(6) < $l(f>(zo), 
so that boundary terms at a, b are neglible, a standard Taylor expansion/inverse Fourier transform 
computation argument near zq [M][p. 60, exercise (ii]) combined with repeated integration by parts 
away from *0 yields 

1(h) ~ h^ +1 )/ 2 e^r^ 

in place of (IB. II) . where j is the order of the first nonvanishing derivative of <f at zo. 

2. Lemma [BJ] might more properly be titled as the Saddlepoint Lemma, as it implicitly involves a 
preliminary deformation of [a, b\ to a curve C passing through a saddlepoint zq, whose contribution 
is then estimated by the method of stationary phase. Note the difference in emphasis between 
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estimate (IB. ID and those of classical stationary phase estimates focused on critical points on [a, b], 
for example, real stationary phase or “Fourier transform” type estimates, or the complex-valued 
version of [PWj . These estimates would typically give a bound of 0(h~ N )e mmh , N arbitrary, 
for 1(h), reflecting absence of critical points in [a, b], but not capturing the full exponential rate of 
decay of (IB.ID . Here, we are focused on the sharp rate of exponential decay, rather than details of 
the asymptotic expansion of the 0(h P 2 ) algebraic factor as in the classical case. Put another way, 
the classical phrasing of the stationary phase lemma as in |PW] corresponds to the case C = [ 6 , a]. 

3. A useful extension, allowing phases 4> = 4>(x, 9(h)) depending analytically on x, 9, 9 monotone 
decreasing with lirrp,^ 0 + 9(h) = 0, is the estimate 


(B.2) 



c/>(x,ip(h)) 

e h 


A(x, h)dx 


, 1 /n (<P(z,):K(h)) + 0(e(h) 2 > 

h ' e h 


(A(z 0 ,0) 



+ 0(h~^ 2 + 9(h ))), 


which follows by first performing the above-described approximate gradient flow for the limiting flux 
4>(-,0) to reduce to the case that zq is a strict minimizer of 5Rp(-, 0 ) on the image of 7 , quadratically 
nondegenerate at t = 0 , then noting that this situation persists under small analytic perturbations, 
by the Implicit Function Theorem and continuity, with the minimizer zj) of (/>(•, h) distance 0(9(h)) 
from the minimizer for <?)(•, 0). We make use of this extension in our computation of (I1.16D above. 
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